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Q_i' Classical vortex solutions in various two-Higgs systems are studied. The systems we 

consider include the standard model with two Higgs doublets, in which case the vortex 
appears as part of a string-like object. The Higgs potentials contain several different 
couplings in general and the spontaneous symmetry breaking involves with two different 
vacuum expectation values. In particular it is shown that the existence of such a solution 
in general requires a specific ratio of the two Higgs vacuum expectation values, i.e. tan (3, 
and some inequalities between different Higgs couplings. This ratio can be determined in 
terms of the couplings in the Higgs potential. The Higgs masses are also computed in 
this case. (l+2)-d solutions are topological so that they are topologically stable and the 
Bogomol'nyi bound is saturated for some couplings. Some comments on the stabilization of 
(l+3)-d solutions are also given. Thus, as long as such a defect can be formed in the early 
universe, stable or not, tan/3 is no longer an independent free parameter in the theory. 

1/93 



* e-mail address: hsla@phys.tamu.edu, hsla@tamphys.bitnet 



1. Introduction 



The idea of the spontaneous symmetry breaking is one of the most important building 
blocks of the electroweak theory, which is one of the greatest achievements in theoretical 
particle physics. Although there is no doubt about the validity of the electroweak theory, 
the final experimental proof however is not yet completed because the wanted Higgs particle 
is still at large. And also there is some possibility that the detail of such spontaneous 
symmetry breaking may be a bit different from the minimal content originally proposed. 
There is growing anticipation, although indirect, the electroweak symmetry breaking may 
be induced by two Higgs doublets, rather than one. 

For example, the recent measurements of the gauge couplings have led us to an- 
ticipation that the minimal supersymmetric Grand Unified Theories (GUTs)0 or super- 
gravity GUTs[Q |4]] with the supersymmetry scale of order 1 TeV or below may be a phe- 
nomenologically plausible unified theory of strong and electroweak inter actions ||. These 
supersymmetric GUTs in general require at least two Higgs multiplets for the electroweak 
symmetry breaking 0] H H- 

Another indication is that the current observed top quark mass bound (mt > 
91GeV) m has already exceeded the upper bound required by the Coleman- Weinberg mech- 
anism of the single Higgs case to radiatively induce the electroweak symmetry breaking||, 
which requires mt < 78GeV|§. Although it is rather sufficient to have tree level breaking 
for particle physics purpose, it will be difficult to have cosmological implications without 
radiative breaking. However in the two-doublet case this upper bound can become suffi- 



ciently high due to the contribution of other scalar field masses ||10|| , so that we can avoid 
such a trouble. 

In exchange of these positive points, having one more Higgs doublet will introduce 
further complication to the theory. Needless to say, first, we have to deal with more 
observable massive scalar particles. Theoretically, it also introduces more free parameters. 
To spontaneously break the SU (2) x U (l)y symmetry down to the U (l) e m each Higgs gets 
its own vacuum expectation value (VEV), say t>i, t>2- These VEVs are phenomenologically 
important but unfortunately they are not determined theoretically except in some no-scale 
models[^TJ. The geometric sum v 2 /2 = v \ + v \ can be determined in terms of the mass of 
the gauge boson, where v denotes the electroweak symmetry breaking scale. This however 
leaves the ratio of the two VEVs, tan/3 = t^/fi, still undetermined. 

If the two-doublet model would turn out to explain the electroweak symmetry break- 
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ing, eventually future experiments will determine tan (3. That will however still leave us a 
question why that ratio should be different from others. This is fairly a common situation. 
We always wonder if nature selects out some particular property among many possible 
choices. Then curiosity drives us to look for some explanation. Pursuing such a question, 
we are often led to a new phenomenon in physics. Thus it is very important to look for 
any argument to constrain the ratio rather theoretically, if possible. 

With such a motivation in mind, in this paper we shall attempt to find any relation 
to constrain tan/3 in two-Higgs systems. The result is indeed positive and we find that 
there is a simple formula to express tan f3 in terms of the couplings of the Higgs potential, 
so far as nature admits certain vacuum defects during the electroweak phase transition. 



Preliminary results were presented injO . 

To demonstrate how it works we shall first work on a simple (l+2)-dimensional toy 
model with U(l) gauge symmetry. Vortex solution in this model takes the role of a nec- 
essary vacuum defect. The generic structure however persists in the two-Higgs-doublet 
standard model, in which (cosmic) string-like solution takes the role.. 

The interest in the application of vortex was initiated in the study of the Ginzburg- 
Landau theory of the superconductivity |13| . Subsequently, in the relativistic field theory 
context, known as the Abelian-Higgs model, the vortex lines are interpreted as string-like 



objects|T4|- Due to the difficulty of solving the nonlinear field equations any exact vortex 
solutions are not yet known in a general case. However, asymptotic solutions can be easily 
found. It is also known that if the gauge coupling constant and the coupling in the Higgs 
potential satisfy a special relation, one can obtain an exact solution [113]. In the Abelian- 
Higgs model one can naively expect that these vortices are stable because of the nontrivial 
topological configuration, i.e. 7Tx(U(l)) 7^ 0. But more careful analysis tells us that this 
is true only if the above two coupling constants satisfy a certain inequality Jl^]. If not, 
solutions in the higher winding sector than one are unstable. 

Such a structure persists even in more realistic models in particle physics. It was 
pointed out that the standard model also admits such a vortex solution (more precisely, 
string-like solution which forms a vortex on a plane perpendicular to the string.) |17j |18j |19] . 



In this case, there is no obvious topological configuration because the hypercharge 
U(l)y is not spontaneously broken, but rather 577(2) x U(l)y is broken to U(l) ern . In 
ref. however, a solution of a pair of SU (2) magnetic monopoles connected by a string is 
derived, (a similar solution can also be obtained in our case. See the section 4.) It was also 
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pointed out that in the standard model case these monopoles can stabilize this system. 
One can also derive string-like solutions in the standard model with one Higgs doublet 
without attaching to monopoles These are unstable solutions except for smdw = 1, 
where 6w is the Weinberg angle.0. 

In this paper we shall extend the results to the cases with more Higgs scalars. If 
there are two Higgs multiplets to start with, we call it a two-Higgs system. As soon as we 
introduce more scalar fields, the potential to induce spontaneous symmetry breaking be- 
comes more complicated. In particular there are more couplings including the interactions 
between different scalar fields. The structure of the phase transition itself may become 
more involved. For example, in a simpler case without gauge couplings and only two scalar 
couplings there are already three different critical points[^T|. If the two VEVs are very 
much different, each scalar can get its VEV one by onei. In the two-doublet model case, 
if so, we can naively anticipate that the cosmic phase transition would occur in two steps. 
Thus perhaps a full investigation of the structure of the fixed points of two-Higgs potential 
may be necessary. 

This paper is organized as follows: In sect. 2 (1 + 2) -dimensional U(l) gauge models 
are considered. One with U(l) x U(l) global symmetry and another without any larger 
accidental symmetry. Then it is shown that these vortices are topologically stable and in 
some case the Bogomol'nyi bound is saturated. In sect. 3 the two-Higgs-doublet standard 
model is investigated. In sect. 4 it is also speculated the possibility of stabilization of 
the electroweak Z-string by attaching monopoles. In sect. 5 the essence of this paper is 
summarized and further issues are discussed. Finally, in the appendix A we show that the 
assumption cyjc^ = V\/v% we made in previous sections is reasonable. 



2. (l+2)-d Models 

As simplified models we shall first consider local U(l) gauge theories with two Higgs 
singlet scalars. In some sense these models can be viewed as a generalization of the Abelian 
Higgs model in the Euclidean two-dimension, but for our purpose which requires sponta- 
neous symmetry breaking, we consider them rather as (1 + 2) -dimensional systems. How- 



1 It has been recently argued that an extra scalar field can stabilize the Z-string if it forms a 



bound state [20 1 . 

2 The t>2 S> vi case was studied in ref. [22] by integrating out the heavy scalar modes. 
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ever solving classical field equations are completely equivalent because we are interested 
in static solutions. 

The Higgs potential we use in fact is motivated by the general two-Higgs potential 
that induces SU(2) x U(1)y —* £7(l) em symmetry breaking, which we shall use in the next 
section and is usually written in terms of two Higgs doublets @ ||. In terms of singlet 
Higgs scalars the potential takes the following form: 

vfa, h) - vi? + ^m 2 - v\f + ^d^i 2 + 1</> 2 | 2 - , 2 ) 2 

\ \ 4 (2-1) 

A 5 
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>l</>2 - VlV 2 



where v 2 = v \ + v\. The significance of the As-term resides in the symmetry structure of 
the Higgs potential. Thus we consider the two cases for A5 = and A5 7^ separately. 

Note that if Ai = A2 = 0, this Higgs system has an accidental SU (2) global symmetry 
if we require that (<pi,4>2) form an SU(2) doublet. This pattern of symmetry breaking, 
SU(2) global x ^(l)iocai — * f (l)giobab is known to lead to semilocal topological defects f23|j. 
If A3 = 0, this becomes simply a decoupled two-scalar system with global £7(1) x £7(1) 
symmetry. There is a vortex solution trivially generalized from the result in ref. [|14| . 

In this paper we shall stick to the general case that A3 7^ and at least one of Ai or A2 
is not zero. Then we shall find that this system reveals a rather interesting result, which 
cannot be obtained otherwise. The key observation is that the spontaneous symmetry 
breaking of Eq. (|2.1] ) leads to a vortex solution, whose existence will introduce an extra 
condition on the Higgs VEVs. Then we can determine them completely, which are related 
only by v 2 = v \ + v\ otherwise. 

Case I: A5 = 

This two-Higgs system has not only the local U(l) symmetry but also (accidental) 
global U(l)i x £7(1)2 symmetry respectively for <pi and </>2-^ If A5 did not vanish, there 
would not be such a global symmetry. In this case this model can also be seen from 
a different viewpoint. One can start from a two-scalar field theory with global £7(1) x 
£7(1) symmetry with a given scalar potential. Then gauge a £7(1) subgroup of the global 
symmetry. Such a viewpoint is common in dynamical symmetry breaking models. 



Such a phenomenon that the potential has a larger global symmetry than the gauge symmetry 



was called accidental symmetry first by Weinberg in [24\ 
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When the two scalar fields get VEVs, not only the local U(l) symmetry is spon- 
taneously broken but also the global U(l) x U(l) symmetry is broken. There are two 
Goldstone bosons due to the global symmetry breaking. One will be absorbed to make the 
gauge boson massive by the Higgs mechanism of the local U(l) breaking. But, as a result, 
there is one Goldstone boson left over. Since the As-term explicitly breaks the global sym- 
metry, we can expect that with As-term it will become a pseudo- Goldstone boson whose 
mass is proportional to A5. This will be considered as Case II. 

Let us now consider the Lagrangian density in (1+2) -dimensional space-time 

C = -\F» V F^ + ±|IV/>i| 2 + i\DM 2 - V(<Pu<P2), (2.2) 

where = d^A u — d u A^ and _D M = <9 M — ieA^. For minimum energy solutions we 
usually solve the minimum energy condition, but in this case it is not so easy to obtain 
the necessary Bogomol'nyi bound. Thus we shall deal with the field equations, then look 
for static finite energy solutions. 

The equations of motion can be obtained from the Lagrangian density as 

D^D^x + (Ai + AsXI^il 2 - vl)<f>! + Xsdhl 2 - vDfc = 0, (2.3a) 

D^fa + (A 2 + AsXI^I 2 - uf)02 + A 3 (|0i| 2 - v\)4>2 = 0, (2.36) 

dT^ = J v = J lv + J 2v , (2.3c) 
J iv = -%ie(<fid v <l>i - fcdvfi) - e 2 A u \^\ 2 , i = 1,2. 

For time-independent solutions we choose Aq = gauge, then the system effectively 
reduces to a two-dimensional one. In this case since we are interested in vortex solutions 
in IR 2 , it is convenient to represent them in the polar coordinates (r, fl) |]14{| such as 

1 = e *»7( r ), cj> 2 = e ine g(r), A = e e -A(r), (2.4) 

r 

where m, n are integers identifying each winding sector. To become desired finite-energy 
defects located at r = these should satisfy the following boundary conditions: 

/(0) = 0, 0(0) = 0, A(0) = 0, 

f —> Vi, g — > V21 A — > const, as r — > 00. 

The constant for the asymptotic value of A will be determined properly later. 
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In the polar coordinates the equations of motion Eq.( [2.3a — c ) can be rewritten as 



--d r (rd r f) + \f(m - eA) 2 + (Ax + A 3 )(/ 2 - v 2 )f + X 3 (g 2 - v 2 2 )f = 0, (2.6a) 



--d r {rd r g) + \g(n - eA) 2 + (A 2 + A 3 )(/ - v 2 )g + A 3 (/ 2 - v 2 )g = 0, (2.66) 
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-d 2 A + ^d r A - e [(m - eA)/ 2 + (n - eA)o 2 ] = 0. (2.6c) 

In general it will be a formidable task to solve these equations exactly, but it is good enough 
to find approximate solutions for large r to show the existence of the vortex solutions. 
Imposing the boundary conditions at large r, Eqs. (|2.6a, fc|) become consistent only if m = n 
and that it fixes the asymptotic value A — >• — as r — > oo. This implies that there is no 
vortex solution of different winding numbers for different Higgs fields. This is in fact an 
anticipated result because the vortex solution we are interested in is due to the spontaneous 
symmetry breaking of the local U(l). With this condition of winding numbers we can solve 
Eq. G2.6dJ for large r to obtain 



A^ - -nJ — ^e r/x + ■■■, 2.7) 
e V 2e 

where A = 1/ev is the characteristic length of the gauge field. Note that the characteristic 
length defines the region over which the field becomes significantly different from the value 
at the location of the defect. 

Now let us determine the characteristic lengths for (pi and <p2 again for large r as 
follows. For simplicity we consider n = 1 case, but the result does not really depend on n. 
Asymptotically we look for solutions of the form 

/ - V! ~ cie" r/Cl , g - v 2 ~ c 2 e" r/ ^ 2 , (2.8) 

where the constant coefficients c\ and c 2 are in principle calculable. The dimensionless part 
of the exact values of c\ and c\ should be at most related to the dimensionless parameters 
of the system or may be pure numerical constants. Thus the essential part of the following 
argument would not be much changed even though we had exact solutions. Also for our 
purpose only the ratio is relevant. Therefore here we shall assume these constants as 
c\ = —v\ and c 2 = — f 2 , which is in a good approximation based on the numerical study 
in the Abelian Higgs model[^5j. Later when we derive the Bogomol'nyi bound, we shall 
also confirm that this is reasonable. (Also see the appendix A.) 
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Then from Eq.( |2.6a"75" ) with Eq.(2TE), in the leading order we obtain 



= Vl e~ r /^ 
= v 2 e~ r/i2 



--j + 2(Ai + A 3 )v? 
si 

-— + 2(A 2 + A 3 )^ 

S2 



+ 2A3^i^ 2 2 e- r /^ + •••. 



+ 2\ 3 v 1 v 2 e 



-r/i 



1 + 



(2.9a) 
(2.96) 



where the ellipses include terms which vanish more rapidly as r — > oo. 

If £i 7^ 4^2: then A3 should vanish. Therefore to have any nontrivial solution for A3 7^ 
we are forced to identify 

Z = ti = 6, 

and that we get the desired result by demanding the vanishing coefficient of e~ r ^ in 
Eqs. Q2.9M ) as 

(2.10) 



tan/3 = — = \/-r-, A3 7^ 0. 
f 1 V A 2 

Thus we have determined the ratio of the two Higgs VEVs in terms of the couplings in 
the Higgs potential. This tells us that although different Higgs field gets different VEVs, 
their characteristic lengths should be the same to form a single defect. Both Higgs should 
reach the true vacuum at the same distance. To do that the two VEVs should satisfy a 
proper relation, which is Eq.( |2.10D . 

Furthermore, together with v, we can completely determine the VEVs as 



t?i = v cos/3 = v\l - — — - r — , V2 = v sin/3 = v 1 



Ai + A 2 

The characteristic lengths £1 > £2 now read 

1 



Ai + A 2 



(2.11) 



£ = 6 = 6 



Ai + A 5 



y/2v V AiA 2 + A 2 A 3 + A3A1 



(2.12) 



Note that, although tan/3 does not depend on A3, it is crucial to have nonvanishing A3 cou- 
pling to obtain such a result. The gauge boson mass is Ma = 1/A = ev after spontaneous 
symmetry breaking. In the next section when a similar structure is applied to a realistic 
model, this gauge field in fact can be identified with a massive neutral gauge boson, e.g. 
Z°. 

In this case one can compute the two Higgs masses to obtain 



mi = 2v 



2 AiA 2 + A 2 A3 + A3A1 2 ^1^2 



Ai + A 2 



m 2 = 2v* 



X1 + X2 



(2.13) 



Note that although the characteristic lengths are the same, the masses of the mass eigen- 
states of the Higgs fields are different. The difference of the two masses is 2X^v 2 . 

The role of the Goldstone boson in this case perhaps need more investigation, which 
we shall leave for future study. 

Case II: A 5 ^ 

The As-term explicitly breaks the global U(l) xU(l) symmetry, so that now the Higgs 
potential has the same symmetry as local U(l) symmetry. Thus there should not be any 
Goldstone boson at all. In fact one can easily see there is another massive scalar field 
whose mass is proportional to A5, signalling that it becomes a Goldstone boson at A5 = 
limit. This is a pseudo- Goldstone boson. 

The As-term in Eq. (|2.1| ) only modifies Eqs.( |2.3a"7£| ). Eq. (|2.3a|) picks up X^((p\(pi — 
V\V<2)<\)2 as Eq. (|2.36j ) does \^,[(\)\(\)2 — viV2)4>i- Thus for A5 7^ we obtain the asymptotic 
equations similar to Eqs.( [2.9a, 6j ) as 

1 



.— + 2(X 1 + X 3 )v 2 1 

si 



+ 2X 3 v lV 2 e- r ^ 2 + X 5 v lV 2 (e" 7 ^ 1 + e" r ^ 2 ) + ■ ■ • , (2.14a) 
=v 2 e~ r/i2 



■ ^+2(A 2 + A 3 K 



1 

si 

+ 2X 3 v 2 v 2 e- r ^ 1 + X 5 vfv 2 (e" r / 51 + e" r ^ 2 ) + • • • . (2.146) 

Since 2A3 + A5 7^ (recall that we assumed all the couplings are nonnegative to obtain 
the minimum at v± and v 2 -), again we have to identify the two characteristic lengths 
£ = £1 = £2 to have nontrivial solutions. Thus the generic structure does not really depend 
on the detail of the Higgs potential so far as A3 > 0. However, there is one more condition 
on A5 we have to impose. From Eqs.( |2.14a, fc| ) for £1 = £2 = £, we obtain 

(Ai - A 5 M 2 = (A 2 - A 5 )4 (2.15) 

Both sides must have the same sign to be consistent. This implies that if A5 is either bigger 
than Ai, A2 or smaller than them, which can be denoted by an inequality 

(A 5 -Ai)(A 5 -A 2 )>0. (2.16) 

From the symmetry's point of view it is more natural to assume the latter case so that 
the As-term would break the global symmetry rather softly. This will also lead to a light 
pseudo-Goldstone boson. 
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The above condition Eq.( 2.l5| ) is again nothing but the constraint on the ratio of the 
two VEVs such that 

Again it does not depend on A3 as before. We can also compute the VEVs and the Higgs 
mass in this case, but the results are equivalent to the doublet case, which will be discussed 
in the next section. 

Remarks on the stability 

Let us first check if the solutions we obtained in this section are topological. This can 
be understood by investigating the topology of the vacuum manifold. In the Abelian Higgs 
model, where there is no accidental symmetry, it is simple to identify the vacuum manifold 
by just looking at the manifold of the equivalent vacuum states of the scalar field. Then 
each topological sector can be identified by measuring the U(l) flux around the vortex. 

In our case if A5 7^ we can also adopt the same philosophy because there is no 
accidental symmetry. If we represent 0i = | <^>i |e z6>1 and §2 = | ^2 1 e* 6 * 2 , then the As-term's 
vacuum condition (p\(p2 — f ii>2 = forces that 9\ = 62- Thus the vacuum manifold is still 
topologically S 1 . For vortex solutions 7fi(S 1 ) = TL so that these are topological solutions. 

If A5 = 0, the topology of the vacuum manifold is a torus, S 1 x S . One may think that 
since tti(S 1 x S 1 ) = 7L © Z, the solutions in this case should be also topological. Without 
the local U(l) symmetry this could be true. However, due to the local symmetry the detail 
is different. Recall that Eqs.( |2.6a, 6p are consistent only if the two winding numbers are 
the same. Thus the topology relevant to the vortex solutions is in fact the diagonal S 1 
of S 1 x S 1 . Although the solutions are topological, but the topology is not completely 
dictated by the Higgs potential. They are dictated by the topology of the local symmetry, 
7Ti (£/(!)) = This is why the two winding numbers must be the same. 

Anyhow, since these solutions are topological, they are topologically stable. However 
as usual it does not seem to be easy to derive the Bogomol'nyi bound for the arbitrary 
couplings. However if the couplings satisfy some relation, we can saturate the Bogomol'nyi 
bound. 

This can be understood more easily if we use the condition £1 = £2, which can be 
rewritten as 

_ 0i _ 02 
V Vi V 2 

Although we have only obtained this condition asymptotically, we shall see this indeed 
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leads to vortex solution. We shall also find out that it is a necessary condition to have 
monopole solution in sect. 4. Thus we believe this will survive as a condition for exact 
vortex solutions. 

We will only check the Case I explicitly, but in all the other cases it can be done 
straightforwardly. With the above condition as we can easily see, the field equations 



Eqs. (|2.3a — 5|) reduce to the field equations of a single Higgs with the VEV v = \/vf + v \ 
and the Higgs potential 

= iA(|</>| 2 - A 

where 

_ A1A2 + A2A3 + A3A1 

A = : : . 

Ai + A 2 

In the single Higgs pointed out in [fl6| , if (3 b = 2A/e 2 = 1, it saturates the 

Bogomol'nyi bound. Thus we have the critical point of the two-Higgs 

2 AiA 2 + A 2 A 3 + A3A1 ro-\a\ 

PB = — 7 — — : = 1- (2.LH) 

e 2 Ai + A 2 

Thus these are stable solutions except for n > 2 if (3b > 1- 



This completes the proof of the existence of vortex solutions in the two-Higgs system 
with the potential Eq.(|2.1|). The consistency condition of such existence led us to be 
able to determine the ratio of the two Higgs VEVs, which otherwise are not completely 
determined. Much work is needed to find exact vortex solutions, but at this moment we 
still get nontrivial physical implication with approximate solutions. 

Note that such a vortex solution in the (1 + 2)-dimensional space is nothing but the 
cylindrically symmetric string-like solution in the (1 + 3)-dimensional space-time. Thus we 
can expect that similar structure should exist in (1 + 3)-dimensional models. This will be 
investigated in the next section. 



3. Two-Higgs-Doublet Standard Model 

One of the most mysterious parts of the electroweak theory lies in the Higgs sector. 
Higgs was introduced to achieve the electroweak symmetry breaking without spoiling the 
consistency of the theory. In this section we shall investigate the structure of a string-like 
defect (so-called "Z-string" ) , which can be formed during the electroweak phase transition, 
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in the two-Higgs-doublet standard model, expecting a similar structure will appear as in 
the £7(1) case considered in the previous section. 

Note that the topology of the local SU{2) x U(l)/U(l) is the same as that of 577(2), 
so iii(SU(2) x £7(1)/£7(1)) = 0. As is well known, the string-like solitonic solutions we get 
would not be topological. This is different from the £7(1) case. However it is known that 
there are nontopological solitons in the standard model[fT7|| [ffS| pi| . 



We shall use the CP invariant two-doublet Higgs potential that induces SU{2) x 
U(1) Y ->• U(l) em symmetry breaking [§ : 

Vfa, 2 ) =|Ai (|^i| 2 - v 2 ) 2 + ±A 2 (|</> 2 | 2 - vlf + ±A 3 (|^i| 2 + \fo\ 2 -vl- vlf 

2 (3.1) 

+ A 4 (\(f>l\ 2 \(i>2\ 2 - l^l 2 ) + A 5 (p\(p 2 - V X V 2 , 

where </>i, 2 are ££7(2) doublets. In this section we shall stick to the general case that 
Xi 7^ for i = 1,2,3 and also assume that all Aj, j = 1,...,5 are nonnegative. This 
potential shows </>i,</> 2 ^ — </>i, ~~ 4>2 discrete symmetry, which is necessary to suppress 
the flavor changing neutral current. Then we shall find that this system reveals a rather 
interesting result, which cannot be obtained otherwise. 

Note that if A4 = = A5, the Higgs potential has a global £7(2) x £7(2) symmetry. 
The symmetry breaking will lead to global £7(1) x £7(1) unbroken so that there will be 
two Goldstone bosons left over. If only A5 is vanishing, due to |0}0 2 | 2 in the A^term the 
global symmetry now becomes £7(2) x £7(1) x £7(1). There still is one Goldstone boson 
after symmetry breaking to global £7(1) x £7(1). 

There however are good reasons to keep A5 7^ 0. First, we don't want Goldstone bosons 
which may not be welcomed phenomeno logically. Secondly, though not directly related, 
in the supersymmetric models A5 is related to the supersymmetry breaking parameter so 
that A5 7^ to have the supersymmetry broken. Note that for A5 7^ there is a global 
£7(2) symmetry, which is locally isomorphic to SU(2) x £7(1). 

Let us consider the bosonic sector of the standard model described by the Lagrangian 
density 

£ = -ItvG^G^ - \F^F^ + \D^\ 2 + \D^ 2 \ 2 - V(<f> u <h), (3.2) 
where = d^B v -d u B^ = d^W* - d u W« + ge abc W*W^ and D„ = d^-ig'^B^- 

a 

ig^W^. Both Higgs' have hypercharge Y = 1. 

Then the equations of motion for the scalar fields are 

=D^D^ X + Ai (|^i| 2 - vl) 0i + A 3 {\<pi\ 2 + \<p 2 \ 2 - vl - vl) 0! 
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+ A 4 (|</>2| 2 </>l - (020l)02^ + A 5 (0201 ~ VlV 2 )(j)2, 

=D»D^ 2 + A 2 (|0 2 | 2 - v 2 ) 2 + A 3 (|0i| 2 + |0 2 | 2 - ^i 2 - v 2 2 ) 2 

+ A 4 ( |0l| 2 02 - (0102)01 ) + A 5 (0j0 2 - ^1^2)01, 



and for the gauge fields we have 



- d^Fau = ju = hv + J2u, 



Jib 



4d v <j>i - (0„&)Vt + k ,2 ^|0,| 2 + y g w^\r a ' 



1,2, 



- d»G% v - ge abc W b »Gc = J- = J*, + J* , 



1*9 



+ Igg'B^r^i + \g 2 W b v 4>\r a T 



(3.3a) 



(3.36) 



(3.4a) 



(3.46) 



1,2, 



For time-independent solutions we choose Bq = = Wq gauge and impose the cylindrical 
symmetry around the string, then the system effectively reduces to a two-dimensional one. 
In this case the string solutions in the (1 + 3)-dimensional spacetime correspond to the 
vortex solutions in IR 2 . When Higgs gets VEV, the false vacuum region forms vacuum 
defects. As usual, we redefine the neutral gauge fields as 



A„ = cos6 w B„ + sin6 w W^ Z„ = sm6 w B„ - cos6 w Wf. 



(3.5) 



where 9w is the Weinberg angle defined by tanOyy = g'/g- We shall also use g = 
\ \/fl f2 + g' 2 f° r convenience. 



Z = e e -Z(r) 
r 



(3.6) 



For vortex solutions it is convenient to represent them in the polar coordinates (r, 6) |1 
such as 

\ , ( 

e im6 fi{r) ) ' 92 \e ind f 2 (r) 

where m, n are integers identifying each "winding" sector ( we shall come back to this 
point later again.). Here we are mainly interested in the case of = = W 2 , but we 
expect there are other solutions similar to the case of ref. ]l 
Then Eqs. (|3.4a, £|) become 



= - -d r (rd r B e ) + \b 6 



r 
9' 



r 



- ~ i( m ~ W B & - 9W$j) fl + (n - \{g'B e - gWl))f 2 ] 
= - -dr{rd r W$) + 



r 
9 



+ ~ i( m ~ W B o - 9Wl)) fl + (n- \{g'B e - gWl))f 2 ] . 

12 



(3.7a) 



(3.76) 



As we can easily see, satisfies a trivial equation so that we can set = 0. Thus from 
the rest of the equations of motion we obtain 



= --drirdrh) + \fi(m - gZ) 2 



+ (Ai + A 3 )(A 2 - vl)h + A 3 (/| - v 2 2 )h + A 5 (/1/2 - <W (n -™ )0 )/ 2 , (3.8a) 



1 



= —d r (rd r f 2 ) + -f 2 (n-gZ) 



2 



+ (A 2 + A 3 )(/ 2 2 - v\)h + A 3 (A 2 - vl)f 2 + A 5 (/i/ 2 - v^e^-^h, (3.86) 
= -d 2 r Z +\z-2g [(m - gZ)f 2 + (n - gZ)f 2 } . (3.8c) 

Note that A4 coupling does not take part in this structure classically. 

To become desired finite-energy defects located at r = the solutions we are looking 
for should satisfy the following boundary conditions: 

/ 1 (0) = 0, / 2 (0) = 0, Z(0) = 0, 

(3.9) 

fi vi, ji —>■ i?2j Z —y const, as r —>■ 00. 

The constant for the asymptotic value of Z will be determined properly later. 

Again we shall look for asymptotic solutions. Imposing the boundary conditions at 
large r, Eqs.( |3.8a, 6|) become consistent only if m = n and that it fixes the asymptotic value 



Z — ► n/g as r — * 00. This implies that there is no vortex solution of different "winding" 
numbers for different Higgs fields. With this condition of winding numbers we can solve 
Eq.( ^.8d) for large r to obtain [p~E 



Z^ --nj — v^e" r/A + ---, 3.10 
9 V 2 9 

where A = 1 /gv is the characteristic length of the gauge field. Note that the characteristic 
length defines the region over which the field becomes significantly different from the value 
at the location of the defect. 

The asymptotic solutions for (pi and </> 2 can be found as follows: For simplicity we 
consider n = 1 case, but the result does not really depend on n. Besides, since these are 
nontopological, it is not really necessary to consider other n sector. Asymptotically we 
look for solutions of the form 

h-vt-cte-^, f 2 -v 2 ^c 2 e- r ^\ (3.11) 
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where Ai and A2 are the characteristic lengths of fa and fa respectively and the constant 
coefficients c\ and c 2 are in principle calculable, thus they are not free parameters. Note 
that we can normalize any dimensionless constants in c, to be the same. Furthermore, for 
our purpose only the ratio is relevant. Therefore these constants can be taken as c\ = —v\ 
and C2 = — v 2 in a good approximation. Even though the exact results differed from these, 
the essential argument of the following is much the same. Then in the leading order we 
obtain 



fie 



v 2 e 



-r/ti 



-r/6 



~ + 2(Xx + X 3 )vl + X 5 vj 

si 

-2 + 2(A 2 + \ 3 )v% + X 5 vf 

S2 



+ (2A 3 + X 5 )v!vle- r/ ^ + ■ ■ ■ = 0, 
+ (2A 3 + X 5 )vfv 2 e- r ^ + ■ ■ ■ = 0, 



(3.12a) 
(3.126) 



where the ellipses include terms which vanish more rapidly as r — > 00. 

Recall that A3 > and A5 > so that 2A3 + A5 7^ 0. Thus to have any vortex solution 
we are forced to identify the two characteristic lengths of the scalar fields such that 

Z = tl= S2- 



We shall find in the next section that to attach monopoles at each side of this string we 
should require that fajv\ = fajv 2 - Thus this is also a necessary condition to stabilize by 
attaching monopoles. 

To be consistent, as in the U(l) case we also obtain an inequality 



(A 5 -Ai)(A 5 -A 2 )>0. 



(3.13) 



However in the two-doublet stabdard model case A5 does not need to be small to be 
natural. Then we get the desired result by demanding the vanishing coefficients of e~ r ^ 
in Eqs. (|3.12a, EJ) as 



tan/3 = — 



X\ — Ar. 

A2 — A? 



A 3 ^0 orA 5 ^0. 



(3.14) 



Thus we have determined the ratio of the two Higgs VEVs in terms of the couplings in 
the Higgs potential. This tells us that although different Higgs field gets different VEVs, 
their characteristic lengths should be the same to form a single defect. Both Higgs should 
reach the true vacuum at the same distance. To do that the two VEVs should satisfy a 
proper relation, which is Eq. (|3.14| ). 
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Furthermore, together with v, we can completely determine the VEVs as 



v v / A2-A5 v . v / A1-A5 , v 

f 1 = — = cosp = — =4 / : — , Vo = — = smp = — = 4 / : — . (3.15) 

y/2 V2\ A! + A 2 -2A5' V2 V% V Ai + A 2 - 2A 5 1 ; 

The characteristic lengths £i,£ 2 , now satisfy 



c = c = c _1 / Ai + A 2 - 2A 5 . 

^ 2 A 1 A 2 + A 2 A 3 + A 3 A 1 -A5(2A3 + A5)' l ' } 

Note that, although tan/3 does not depend on A3, it is crucial to have nonvanishing A3 
or A5 coupling to obtain such a result. The gauge boson mass is Mz = 1/A = gv after 
spontaneous symmetry breaking. 

In this two-Higgs-doublet model there are five physical Higgs bosons: H , A , if , h°. 
A is a CP-odd neutral scalar, while if , /i° are CP-even scalars. h° denotes the lightest 
Higgs. Using Eq.( |3.15| ), we can compute the masses of all these physical Higgs bosons in 
terms of the couplings in the Higgs potential and v, where v — 247 GeV. A5 is related to 
M^o and MH ,h can be determined in terms of Ai, A 2 , A3, A5 and v. Thus we only have 
five free parameters, if such an electroweak Z-string exists. 

In particular the neutral Higgs masses become impressively simple: 



mffo = v 2 



AiA 2 -A 5 

A 3 + 



m/jO = v 



Ai + A 2 — 2A5 
2 2AiA 2 — As(Ai + A 2 ) 
2(Ai + A 2 - 2A 5 ) 



(3.17) 



Note that m h o does not depend on A3. 

The appearance of integers in the solutions, which we still call "winding" number, is 
rather intriguing because there is no explicit U(l) symmetry to be broken which should 
determine the necessary topological sector. If our vortex solutions are nontopological as in 



ref. [2(J, there should not be such a parameter. This however can be explained as follows: 
If we regard W} = = W? as gauge fixing conditions, then effectively we can view the 
symmetry of the system as U (1) x U (1)y ■ When we twist this symmetry to obtain £7(1) em , 
the remaining twisted U(l)~ is spontaneously broken to lead to the winding sector. 

Since SU(2) is a simple group, U(l) x U(1)y is not an invariant subgroup of SU{2) x 
U{1)y- Furthermore, tci (SU(2) x U(1)y /U(l) em ) = implies that these winding sectors 
would not provide any topological stability. In other words, they must be gauge equivalent 
to n — 1 solution via deformable gauge transformation. 
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Even for n = 1 solution it is most likely that this solution would not saturate the 
Bogomol'nyi bound. Although it obviously is a finite energy solution, it does not seem to 
be a classically stable solution. It is argued that there is a case of quantum stabilization 



of a classically unstable solution J27|. It however cannot be applied in this case unless one 
of A4 or A5 vanishes because there is no tree level Goldstone boson. If A5 = 0, then there 
is one Goldstone boson. But the argument still does not apply because 7ri(G , g au ge/-ffg aU ge) 
is still trivial. 

Nevertheless, there is one more possibility to stabilize such an electroweak Z-string, 
which will be presented in the next section. 



4. Monopole-String-Antimonopole 

In the one-doublet standard model it is claimed that such a string solution can be sta- 
bilized by attaching monopoles at each end and keeping them sufficiently far apart [ff7|. The 
rationale behind the stabilization is that these monopoles are genuine SU(2) monopoles 
so that they can be topologically stable. In other words, the string can be constructed as 
a string connecting two Wu-Yang type monopoles. 

In fact we can apply the same argument here. As we shall see soon, the two-doublet 
model also admits such a monopole solution. The consistency condition for the existence 
is again <pi/vi = ^2/^2, which is an equivalent condition to £1 = £2 in the string case. 
Thus it seems to us that this condition is not just accidental but may have more profound 
significance in the spontaneous symmetry breaking of the two-Higgs systems. 

In the static case the energy density of the system can be easily derived from the 



Lagrangian density Eq. (|3.2|) . Then the minimum energy conditions are 



Dfj,(f>i = 0, D^ 2 



and 



Eq. (O) implies that 



V(<j> u <h 



|0i| 2 = «?, Ifol 2 = v%, 4>\<f> 2 = n r 2 . 



(4.1) 

(4.2) 
(4.3) 



where the third condition does not follow from the two other conditions. 

From Eq. (|4.1|) by multiplying </>Jr a and subtracting its hermitean conjugate we obtain 



1,2. 



(4.4) 
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These two simultaneous equations need to be consistent to have solutions. Comparing 
them, we are required to identify 

^ = (4.5) 
Vi v 2 

Thus the two equations become identical. Then using the Nambu's method we can solve 
Eq.(fOj) to obtain and _B M in terms of cf>i and fa. 

For a monopole solution with a singularity along the negative ^-direction, the Higgs 
fields would be 

<k = Vi[ • xl iip ), r^O, 1 = 1,2, (4.6) 
where the singularity is identified by the ill-defined phase at 9 = n. Thus we have 

4r a ^ = v^y. (4.7) 
Since we expect that behaves like a monopole with a Dirac string, we use 

9'B a = -e^^— v (4.8) 
r[r + z) 

where we have identified the spatial indices and the gauge group indices as a, 6, ... = 1, 2, 3 
and In the spherical coordinates we can easily see that only 

(^-component is nonvanishing. Bq = by the gauge choice. 
Then plugging into Eq. (|4.4|) we can solve for to get 

gW b = _ e a bc *_ (49) 

which is nothing but the SU{2) monopole solution[EM|. Note that (j)\r a (j)i takes the role of 
the adjoint Higgs field in the usual non-Abelian Higgs model. Since denotes a genuine 
SU(2) monopole, the string carries only £7(1) returning flux but does not carry SU(2) 
returning flux. 

Now we can patch a pair of monopole and antimonopole, which leads to 



cos ^0 
sin±6e^ 



(4.10) 



where cosG = cos6*i — cos#2 + 1 and 6>i (62) is measured from the position of the monopole 
(antimonopole). If the distance between the monopole and the antimonopole is sufficiently 
far apart, the solution we derived can be used for each side so that this can be viewed as 
a monopole-string-antimonopole system. One can easily obtain the relevant solutions by 



generalizing Eqs.( [4.8| ) ( fl.9| ) in this case. Thus by the same reasoning as in |T7[], this would 
be stable. 
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5. Discussion 

We have shown that two-Higgs systems in general admit vortex solutions, which re- 
quires a specific ratio of the two VEVs. This can be understood more easily if we use the 
condition 

J= ti = 

Recall that this condition not only shows up as a consistency condition of the vortex 

solutions but also shows up as that of the monopole solution considered in sect. 4. 

Then as we can easily see, the field equations reduce to the field equations of a single 

Higgs with the VEV v and the Higgs potential V{<p) = \(\4>\ 2 - f 2 /2) 2 , where = v<p/y/2 

(in the 2-d case replace v 2 /2 by v 2 ) and 

_ A1A2 + A2A3 + A3A1 — As(2A3 + A5) 
Ai + A2 — 2A5 

All such results can be obtained if 

tan^^ = ,/5H^, A 3 ^0. 
ui V A 2 - A 5 

These results also apply to the cases of A 5 = 0. 

Although there are characteristic differences, in fact we have observed that so far as 
vortex solution is concerned, such a structure is quite generic for both cases: the (1+2)- 
dimensional £7(1) system and the (l+3)-dimensional SU(2) x £7(1). Two big characteristic 
differences are that the former is topological and saturates Bogomol'nyi bound (in some 
case), whilst the latter is nontopological but can be stabilized by attaching monopoles. 

In this paper we have dealt with the tree level potential. It is usually believed that 
tree level solution reappears in loop-corrected effective potential, although there might be 
quantitative differences. This is in fact relevant to study the porperty of vacuum defects 
in the context of the cosmological phase transition. 

For the unstable solutions we need further investigation to find out what kind of 
cosmological trace they could leave. But one certain thing is that as soon as we find out 
how nature selects out tan (3 (if nature prefers the two-doublet model), we shall understand 
its mystery if tan/3 turns out to meet our claim. The ultimate proof of the existence of 
such vacuum defects should be determined by experiments or by observations. 

It will be also important and interesting to find out if there is any other reason nature 
prefers a specific tan (3 theoretically. This will be still a question to be answered after we 
determine experimentally. In this sense, we hope this work can provide a clue to future 
investigation. 
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Appendix A. 

In this appendix we shall show explicitly why C1/C2 = vi/vz is a reasonable assumption 
in the case I of the section 1. In other cases it can also be generalized straightforwardly 
so that we leave them as simple exercises. The main rationale behind our assumtion is 
that these constants are not usually determined by solving nonlinear differential equations 
asymptotically so that one should use them as input parameters. Of course they will be 
determined by solving equations exactly, however we are not yet able to solve the equations 
in question exactly. As we shall see, these will be ill-defined if we try to determine otherwise 
and the assumption we made is the only reasonable consistent choice. 

In eq. (|2.8|) we can conveniently redefine q = c^, where q are dimensionless. Also 
let c = Ci/c2, then the condition for = £2 reads as 



Suppose c is not an input parameter but to be determined by this equation for any tan (3 
to nullify our claim that this equation relates tan/3 and other couplings, then one should 
be able to determine c for given tan (3. In this case eq. ( |A.1[ ) becomes a quadratic equation 
for c as 

A3C 2 + (tan 2 (3{ A 2 + A 3 ) - (Ai + A 3 )) c - A 3 tan 2 (3 = (A.2) 
and one can obtain 
1 



c± = 



2A, 



- tan 2 /3(A 2 + A 3 ) + (Ai + A 3 ) ± \J (tan 2 /?(A 2 + A 3 ) - (A x + A 3 )) + 4A 2 tan 2 /? 

(A.3) 

First, note that c_ < so that one of the Higgs fields approaches to the true vacuum 
from the wrong direction. Furthermore, c _ — > for any A's as tan/3 — > 0, so c_ leads to 
ill-defined coefficients. Thus we are left with c+ > 0. 

Second, if A3 — ► 0, then any c is ill-defined unless the coefficient of the linear term in 
eq.( [A.2| ) vanishes. This vanishing condition is nothing but eq.( [2.10| ) as A3 — > 0. Otherwise, 
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c either increases indefinitely or approaches to indefinitely. This is unreasonable because 
for vortex solutions with the same characteristic length we should expect c\ and C2 are of 
the similar order for any A3. Without loss of generality, say, ci <C C2 ~ lj then the notion 
of characteristic length in the asymptotic formula for (pi field we assumed fails to make 
sense. 

Therefore, eq. ( |A.2| ) cannot be solved for plausible c consistently for tan (3 as an input 
parameter. This leads us to a conclusion that we should treat c as an input parameter and 
the assumption we make in the previous sections is reasonable. Thus c = 1, which leads 
to eq.( |2.10| ) consistently, is a well-defined reasonable assumption to present our argument. 
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